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Abstract. - We derive a closed equation for the shape of the free surface of a magnetic fluid 
subject to an external magnetic field. The equation is strongly non-local due to the long range 
character of the magnetic interaction. We develop a systematic multiple scale perturbation 
expansion in which the non-locality is reduced to the occurrence of the Hilbert transform of the 
surface profile. The resulting third order amplitude equation describing the slow modulation 
of the basic pattern is shown to be purely local. 



Introduction. - The emergence of spatio-temporal order in distributed systems can often 
be theoretically analyzed in terms of amplitude equations ^E]. Generically these amplitude 
equations are non-linear partial differential equations for the slow time and space variations 
of an envelope function of unstable modes. As such they reflect the local character of the 
underlying dynamics. 

In several interesting cases, however, the dynamics of the system is substantially influenced 
by non-local interactions. These may arise, e.g., from a mean flow in convection problems 
EJ 03 El) from electric and magnetic fields in solids 0EI or electrochemical systems [§], 
from long range elastic interactions |10j . or from global couplings in systems with chemical 
reactions |lllll2llT3] . 

Since the type and stability of the emerging patterns is usually modified if non-local in- 
teractions are present the corresponding amplitude equations are expected to exhibit some 
degree of non-locality as well. From a phenomenological point of view one may therefore 
be tempted to simply add to the standard form of an amplitude equation non-local terms 
complying with the relevant symmetries of the system |14U15| . A systematic derivation of the 
amplitude equation from the basic non-local dynamics has been accomplished in a few cases 

only. 0[miMMni] 

In the present letter we investigate the formation of static patterns on the free surface of 
a magnetic fluid subject to an external magnetic field. In this process the non-local character 
of the magnetic interaction is of vital importance. We establish a closed, strongly non-local 
equation for the free surface of the magnetic fluid and systematically derive an amplitude 
equation for the surface deflection in the vicinity of the critical field strength. Although the 
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corresponding linear operator is non-local the lowest order amplitude equation is found to be 
local and to reproduce previous results from a perturbative analysis of the free energy |18U19| . 

Magnetic fluids are suspensions of ferromagnetic nanoparticles in suitable carrier liquids 
behaving as super-paramagnetic Newtonian liquids If the strength of a magnetic field 

perpendicular to the flat free surface of a magnetic fluid exceeds a threshold value the surface 
is known to become unstable due to the normal field or Rosensweig instability JJj and a 
hexagonal or square pattern of fluid peaks develops |17II18| . 

In order to keep the analysis simple it is convenient to consider a situation with just one 
unstable mode analogous to the roll solution in hydrodynamic systems. As was known experi- 
mentally for some time [30] and has been clarified theoretically recently [2J the corresponding 
"rigde" pattern on a magnetic fluid surface may be induced by an oblique magnetic field. In 
this case the field component tangential to the flat surface suppresses surface deflections in 
this direction. 

Basic equations. - We consider an incompressible magnetic fluid of infinite depth, density 
p and constant permeability fi = /iofi r in a magnetic field Ho which, in the absence of any 
magnetically permeable material, is homogeneous and of the form Ho = Hze z + Hxe x - The 
free surface between the fluid and the magnetically impermeable air above is described by 
z = ((x,y) (see fig. QJ. The surface tension is denoted by a and gravity acts parallel to the 
z-axis, g = — ge z . We assume the horizontal magnetic field component Hxe x to be strong 
enough to suppress surface deflections in ^-direction rendering the problem effectively two 
dimensional. 

The static surface profile Q(y) is then determined by the pressure equilibrium at the surface 



d y ( /i (/x r -l) 2 /io(Mr-l) 
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Here H| and H n \^ denote the magnetic field and its normal component, respectively, in the 
fluid at z = C(y)- The air pressure p on the right hand-side is assumed to be constant. The 
first two terms of eq. Q are the standard expressions describing the influence of gravity and 
surface tension. The remaining terms characterize the impact of the magnetic field on the 
surface profile C(j/)- 

The magnetic field H(y, z) has to fulfill the magneto-static Maxwell equations 

V • B = and V x H = 0, (2) 

where the magnetic induction B(y, z) = ^o/i r H(y, z) is a linear function of the magnetic field. 
The field equations are supplemented by the boundary conditions 

(B (a) -B (f) )| c n =0 (3) 

(H( a >-H< f ')| ( xn =0, (4) 

which describe the feedback of the surface profile on the magnetic field. Here 

_ (o,-fyC.i) 



(5) 

is the normal vector on the free surface pointing outwards and the upper indices ( a ) and ( f ) 
refer to the air above and the fluid below the interface, respectively. The conditions 

H(»,z)-( „ HzGz ^ Z ^°° (6) 

v ; 1 H z /fi r e z if z -> -oo w 
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Fig. 1 Fig. 2 



Fig. 1 - Static surface profile ((y) and magnetic field lines for a magnetic fluid with /j, t = 2 in an 
external magnetic field corresponding to e = 0.0725. The critical wavelength A c = 2n/k c is typically 
of the order of 1 cm. 

Fig. 2 - Stability diagram for a magnetic fluid with = 2 in an overcritical magnetic field, e is 
related to the magnetic field strength by (JSJ, S denotes the deviation from the critical wavenumber 
k c = 1. N is the neutral curve of the linear instability, E denotes the Eckhaus boundary following 
from the amplitude equation I29H . Also shown is the wavenumber fc m i n minimizing the free energy as 
determined in |19|. 



ensure that the magnetic field asymptotically approaches its externally prescribed values. 

It is convenient to rescale lenght by the capillary length fc^ 1 = s/cr/pg, pressure by y/J 
and magnetic field by the critical field strength of the Rosensweig instability 



y MO (Mr - I) 2 

Moreover we introduce the supercriticality parameter 

e = Hi - 1 

and the parameter 

V 



Mr - 1 



(7) 



(8) 



(9) 



Mr + 1 

characterizing the magnetic susceptibility of the fluid. Denoting the tangential component of 
the magnetic field by Ht eq.JQ) then assumes the dimensionless form 



c- 



(c\C) 2 ] 3 / 2 v 1 
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Mr tt2 
H n 



1 



n 



(10) 



where the value of p was determined from the reference state £ = 0. 

The dependance of the magnetic field on the surface deflection via the boundary conditions 
iJSJ and Q is rather implicit. To make the non-locality of the problem more explicit we take 
advantage to the linearity of the magnetic field problem and employ a formal Greens function 
solution for the magnetic field. To this end we write the field in the form 



H(y,2) 



Mr - 1 
2Mr 



Hz 



dy' , i \ (y-y')e y + (z-C(y'))e z l + Mr 
2n Q{y) {y-y'Y + (z-ay')) 2 ^ 
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where the so far unknown magnetic source density q(y) is nonzero at the surface z = £(y) 
only. The field already fulfills the Maxwell equations as well as the boundary condition 
0. The remaining boundary condition Q gives rise to an equation for q(y). To derive this 
equation we observe that the normal component of the field (jlll) is discontinuous at the 
interface, 



Mr 



1 



(12) 



Moreover, interpreting the integral in (|llfl for z = ((y) as Cauchy principal value we find for 
the field directly at the surface H n (yX(y)) = (Hn + Hn^)/2. Together with Q this gives 
rise to 



and using this relation in eq. I|llfl we obtain for q(y) the linear integral equation 

q=l- r)(d y QSq + V Tq. 
The operators S and T act on bounded functions f(y) and are defined by 



+oo 



(Sf)(y) = 



y-y 



dy' 

V J^y'Y + (((y) - Civ 1 )) 2 



(13) 



(14) 



(15) 



and 



(Tf)(y) = i d 4- 



C(y) - C(y') 



* (y-y') 2 + (ay)-ay')) 2f{y,) ' 



(16) 



where j- denotes the Cauchy principal value of the integral. Note that due to the constant 
term on the rhs of l|14|l also the asymptotic conditions (jfJJl are fulfilled. 

With the help of (|11|) we can express the magnetic field contributions in (|10|l in terms of 
the source density q(y) and obtain the pressure balance in the form 



o = C 



d 2 X 



1 + e 



[i + (9,0 2 ] 3 / 2 v i + (^C) 5 



fr[i + (a„C) 2 ]Sg + (a„C)g) , g 2 
(1 + 17)7/(1-17) 77 



(17) 



Eas. l|17|l and (|14fl are a closed system of integro-differential equations for the stationary 
surface profile £(y) and the related source density q(y). 

While the kernels of the operators S and T are nonlinear functions of the surface deflection 
C(y) eq. (|14|) is a linear equation for the source density q(y). For given £(y) it can hence be 
solved formally using the Neumann series 



q = l + J2r 1 n [T-(d v OS} n l. 



(18) 



Replacing q(y) in eq. (jT7|l by this series we obtain a closed nonlinear integro-differential equa- 
tion for the free surface profile C(y). This equation contains products of an infinite number of 
integral operators and is per se of limited use only. Nevertheless it forms a suitable starting 
point for a systematic multiple scale perturbation analysis of the weakly non-linear regime. 
On the other hand this expansion can also be performed directly on the system 11170 and (|14|) 
which is what we do in the following. 
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Weakly nonlinear analysis. We consider slightly supercritical fields < e <C 1 and 
assume that both the surface deflection £(y) and the source density q(y) can be expanded in 
powers of e 1 / 2 : 

C = £ 1/2 Ci + e C 2 + e 3/2 C 3 + ... (19) 

q=l + e 1/2 q 1 + eq 2 + e 3/2 q 3 + .... (20) 

For the resulting expansions of the linear integral operators S and T it is convenient to start 
with the series representations 

Knr N _ >fi 7V ffi W - CG/)) 2 "/Q/)] 



n=0 



and 



(T/)(y) -^ o (2n+l)!/ n (y> - y) ' (22) 

arising from an expansion of the denominators in 115(1 and ((161) and appropriate partial inte- 
gration. The central building block in these expansions is the Hilbert transform defined for 
bounded functions h(y) by 

— oo 

Using CJjJ 

in eas. l(21() and ((22(1 gives rise to expansions of S and T in powers of e 1 / 2 which 
together with 1(19(1 and 1(20(1 are used in ((17(1 and ((14(1 . Matching powers in e then yields a 
hierarchy of linear equations for the expansion coefficients of £ and q in the usual way. 
To order 0(e 1 / 2 ) we find 

CCi -Ci+2a y HCi-9 2 Ci -0. (24) 

The linear operator £ is non-local due to the occurence of the Hilbert transform originating 
from the non-local magnetic interaction. With the help of TLexp{iky) = sgn(fc) i exp(iky) the 
solution of this equation is easily found to be 

Ci(y) = Ae iy + A*e- iy = Ae ly + c.c. , (25) 

with a complex valued amplitude A. As it should the first order of the perturbation expansion 
hence reproduces the results of the linear stability analysis, in particular k c = 1. Note also 
that under the standard scalar product the operator C is self-adjoint, £ T = C. 

For e > not only the critical mode but a whole narrow band of modes with fc-values 
around k c are unstable. As is well known this gives rise to long wavelength modulations of the 
basic pattern [Q . To take account of this possibility we allow for a slow spatial dependence of 
the amplitude, A(Y), formalized by the substitution d v — > (d y + e^l^dy). Using the properties 
of the Fourier transform of A(Y) one can also show TtA(Y) exp(iky) = A(Y)7i exp(ifcy). 

To order O(e) we then find 



£( 2 = 2rjA 2 e 2ty + c.c. , 



(26) 
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with the solution 

( 2 {y) = 2 V A 2 e 2ly +c.c. (27) 

Here the amplitude of the solution of the homogeneous equation was set equal to zero in 
accordance with the general requirement that all Cn with n > 1 should be orthogonal to £i . 
Eventually, to order 0(e 3 ^ 2 ) we get the inhomogeneous equation 



£( 3 = 2Ae iv 



d Y Ae l 



§-«*" 



\A\ 2 Ae iy + \ir\AdYAe 



2iy _ 



A 



3 g3iy 



C.C. 



(28) 



Since C is singular (cf. eq. iJUJl) the right hand side of this equation must be orthogonal to 
the zero eigenfunctions of & which are e ±ly . Returning to unsealed variables we hence find 
the amplitude equation 



= eA+-d 2 A- 



5 A 2 
4rr 

4 ' 



\A\ 2 A. 



(29) 



Despite the non-local character of the linear operator C this equation is local. The terms 
without derivatives are identical with those obtained from a perturbative analysis of the free 
energy of the static surface deflection 18, 19 . In particular we reproduce the result that for 
■q > v5/4 or equivalently fi r > 3.5353 the cubic term becomes positive and higher order 
non-linear terms are needed to saturate the instability |23II24| . 

After imposing the solvability condition l|29|l we can integrate l|28[) to get 



( 3 ( y ) = 4i v Ad Y Ae 2ly 



Ajf 



A 3 e 3iy 



(30) 



Eqs. I|25|) . H27fl . and (|30|l with dyA — have been used to make the plot of £(y) shown in 

Fig.m 

Building on the dispersion relation for surface waves on ferrofluids (25U26| also a slow time 
dependence of the amplitude A with the corresponding term on the left hand side of 1)29(1 can 
be included. This then allows e.g. to investigate the stability of solutions with slow spatial 
modulations of the form 

A(y) = Ae iS y. (31) 
A straightforward analysis £Q yields as boundary for the associated Eckhaus instability 

5 2 E (e) = 2e/3, (32) 

which is shown in Fig. [21 

Conclusion. To determine the equilibrium shape of the free surface of a magnetic 
fluid in an external magnetic field is a highly non-local problem. Eliminating the magnetic 
field problem by introducing a suitable magnetic charge density on the free surface we have 
derived a closed set of two non-local equations for the determination of the surface profile. In 
a systematic multiple scale perturbation expansion of these equations the non-locality results 
in a non-local linear operator C involving the Hilbcrt transform. On the other hand, the 
third order amplitude equation for the surface deflection was shown to be local. An extension 
of the perturbation expansion to other patterns as hexagons and squares or higher orders is 
straightforward though tedious. 
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